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SUMMARY 

A previous report gave calcxLLations for the pressure distri- 
bution over thin obliqiie airfoils at supersonic speed. The present 
report extends the calculations to subsonic speeds. 

It is found that the flows again can be obtained by the super- 
position of elementary conical flow fields. In the case of the 
swept— back wing the pressure distributions remain qualitatively 
similar’ at subsonic and supersonic speeds. Thus a distribution 
similar to the Ackeret type of distribution appears on the root 
sections of the swept— back wing at M = 0. The resulting positive 
pressure drag on the root section is balanced by negative drags on 
outboard sections. 


INTRODUCTION 

So far as is known, attempts to extend airfoil pressure- 
distribution calculations to three-dimensional flow have been 
confined to cases of thin lifting surfaces. It has generally been 
assumed that the component of the pressure distribution arising from 
the thickness of the airfoil will be but little affected by the 
finite span, or aspect ratio, of the wing. This supposition is 
borne out by the known incompressible— flow solutions for flat 
ellipsoids. These solutions show that the usual variations of 
aspect ratio produce small effects. 

Compressible— flow theory shows however that the effects of 
plan form become more pronounced at hi^er speeds. The theory 
Indicates a. progressive reduction of the equivalent aspect ratio as 
the Mach number approaches 1.0. Hence at these speeds the three- 
dimensional character of the flow can no longer be neglected. Of 
particular interest are the deviations from two-dimensional flow near 
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the root sections of a swept— hack wing, since the adverse effects of 
compressihility may arise first in this region. 

In the present report three-dimensional flows are obtained from 
a distribution of "pressure sources" in the chord plane of the airfoil. 
The shapes thus obtained are symmetrical airfoils at zero lift. 

The calculations are simplified by considering airfoils composed of 
conical or cylindrical surfaces In these cases the sources can be 
arrar.ged into lines of uniform strength following the generators of 
the surface. The relation between the strengths of the line sources 
and the shape of the airfoil is the same as in reference 1, that is, 
each line source produces a. deflection of the streamlines crossing 
over the source. The pressure field of the line source again can be 
represented by systems of straight rays of equal pressure (Isobars) 
radiating from the ends of the line soiu:*ce. 

In general, the present development follows closely that of 
reference 1 and the reader should consult that report for addi- 
tional details of the method. The solutions are given explicitly 
for M = 0 but are extended to other Mach numbers by the well-known 
Prand tl transf ormat i on . 


THE OBLIQUE LDIE SOURCE 

It is well known that an individual velocity component of a 
potential flovr will satisfy the same differential equation as the 
potential. In the approximations of the th in-airfoil theory the 
pressure depends only on the individual component u, that is. 



( 1 ) 


while the slope of the surface depends only on the individual 
component w, that is 


dz _ w 
dx V 


( 2 ) 


(Se© appendix for symbols.) Hence in the thin— airfoil theory it is 
often more convenient to deal directly with the velocities u and 
w as solutions of Laplace's equation than to derive these components 
from a velocity potential 

Since u is proportional to the pressure, a solution of Laplace's 
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equation can represent directly the pressure distrihution, hence 
the terra "pressure potential." In this terminology, the fundamental 
solution 



1 


V 


P P P 


(3) 


represents a. imit point source of pressure rather than a point 
source of fluid. 

To get the effect of a row of sources, or a line source, along 
the x-aocis between the points a and h it is necessary to integrate 
equation (3): 

u = / = 3inh~^ — sinh~^ (4) 

Ja v^(x-|)^+y^+z^ Vy^+s^ J y^+z^ 


The pressure field of the finite lino source thus consists of the 
sum of tvro conical pressui'e fields radiating from the ends of the 
line source. (See fig. 1.) In the supersonic case (reference l) 
the radial isobars forming the conical field were confined to the 
downstream Mach cone. Here however the isobars extend over the 
whole space. ^ 

If the direction of flight is along the axis of the sotirce 
(x— axis), the flow will satisfy the boundary condition for a body 
of revolution. However, if the line source is turned out to a 
position oblique to the stream, the boundary shape will be distorted 
and, if the angle of obliquity is large enough to place the line 
source well outside the diameter of the original body, the figure 


The conical pressure field for either the subsonic or the super- 
sonic line source may be obtained directly from the general solU'=- 
tions of Laplace’s equations of zero degree in x,y,z given by 
¥.F. Donkin. (See reference 2, page 357-) The general solution is 

u = f( 

^ x+V x^+y^+z^ y 
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formed will be an oblique wedge. The nose angle of the wedge is 
formed where the streamlines of the main flow cross the line source. 

At supersonic speeds the expiression for the oblique line source 
was obtained by applying an equivalent of the Lorentz transformation 
for which the wave equation is invariant. The equivalent trans- 
formation for Laplace's eqmtion is a rotation of the axes, given by 


X* = X + my 
y' = y - mx 
z' = (v/l -f m^ z 

where m is the slope of the new axes relative to the old. (Note 
that a change of scale is admitted for convenience . ) The geometrj'’ 
of the pressure field relative to the line source is not altered in 
any way by this rotation an'i the isobars behave as thou^ they were 
rigidly attached to the ends of the source. For a line source with 
one end at the origin we have 


u = ainh~^ — (5) 

yo' 

This field is illustrated in figure 2 for the plane z = 0. As 
m — > 00 the X— and y— axes interchange and there is obtained 

u = sinh“^ ^ (6) 

V x2+z2 


for a line source along y. 


See footnote 1, pg. 3- 

The solution corresponding to the subsonic line source is 

y^-lz 


u = — E.P. log 


= sinh' 


x+v' x‘^+y‘=+z‘^ ,/y2+z^ 

while the field for the supersonic source is given by 


u = - E.P. log =z^==== = coah ^ 

xW x“^— y — z 


J 


.r^+z‘ 
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The vertical velocity w near z = 0, which determines the 
shape of the boundary, may he found by integrating u with respect 
to X and then differentiating the resulting velocity potential 
with respect to z 


Evaluation of this integral for the overlapping fields from two 
ends of a line source gives 


over the area of the x, y plane behind the line source . 

The figure formed by the streamlines crossing a line source is 
thus a. wedge-shaped body having an oblique leading edge and extending 
indefinitely downstream. It is evident from equation (3) that the 
infinitely wide wedge cannot be treated in subsonic flow, since it 
creates an infinite presstire disturbance at all points. 

The slope of the wedge surface a.way from the chord plane is 
given by 


With this relation and equation (8) the pressiire coefficient near 
the plane z = 0 may be ea^jressed in terms of the slope 



(7) 


vr = +2Jt 


( 8 ) 


m 


zXP _ 1 m dz 
q 3C dx 


( 



( 10 ) 


where |y’j indicates the absolute magnitude of y* . Following the 
thin-airfoil theory, the pressiur-e over the chord plane (z — 5>0) is 
taken as the pressure over the actual airfoil surface. 
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AIBFOIIS BOUNDED BY PLANE SUEPACES 

It was seen that the effect of a line source in the pressure 
field is to cause a. deflection of the streamlines crossing the source. 
The deflection thus produced is equal and opposite at points above 
and below tho chord plan(^ so that the source spreads the streamlines 
apart. If tho soiurce is followed by a sink of equal strength, an 
equal opposite deflection of the streamlines will occur as they 
cross over the sink. The figure formed by the streamlines near the 
plane z = 0 will thus be a plate of uniform thickness with a 
beveled loading odge.^ (See fig. 3) 

Tho pressure distribution over su.ch a beveled edge may be 
obtained very simply by superimposing the pressures laid off on 
radial isobars originating from tho four comers of the bevel. 

Figure 3 illustrates this process for a bevel having a. square plan 
form. Only isobars from one tip are shown because of the symmetry 
of the figure. 

In figure 3 the line source and the line sink are parallel to 
tho y-axis, hence 


u = sinh ^ ■ / •- - - - sinh ^ - ainh + s-^nh ^ 7^"^ (H) 

|x+l| |x+ll jx-ll [x-l] 

It can be soon that if the aspect ratio of the figure is increased 
to a large va,luc tho ends of the line sources will be separated by 
a great distance and the isobars in intermediate regions will 
approach parallel straight lines, hence the flow field approaches 
a cylindrical or two-dimensional form. At the same time the 
arguments y±l/|x±lj in equation (11) become y+q/xil and q 
takes on very largo values so that 


s 


i„k-i /±a_ 

|x±l| 


— 5> ± log 2 


y±q 

|x±il 


2According to the thin-airfoil theory the thickness of the figure 
ends abruptly at the ends of the source lines . A more exact 
consideration would bo expected to show some rounding at the tips 
of the wedge as indicated in figure 3* 
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and equation (11) is found to approach the Legendre function Qo^ 
that is 


u = 2 log 


x-1 

x+1 


= - 4 Qo (x) 


( 12 ) 


(SoQ reference PS- 110.) 

This expression when comhined with equation (8) agrees with the two- 
dimensional potential function for the wedge, that is. 


- (u-iw) = 4 Qo (x)±2jti Pq (x) (13) 


(See fig. 4.) 

The isobars at ri^t angles to the axis of the line source are 
lines of zero pressure, hence the rays originating at the tip of a 
rectangular wing contribute nothing to the pressure distribution at 
this tip. The whole pressure distribution at one tip is thus obtained 
by considering only those isobars radiating from the opposite tip. 

It is evident that in the case of a long narrow rectangular wing 
the pressures at either tip will be approxima.tely one-half the 
pressm’es over the middle portion of the wing. 

In case the wing is oblique the tip sections will no longer 
be at right angles to the axes of the sourco lines and the rays 
originating from the adjacent ends of the source lines will contri- 
bute to the pressure over the tip. It can be shown that this ccanpo— 
nent of the tip pressure distribution is similar in form to the 
Ackeret typo of distribution, that is, the pressure at any point 
of the surface is proportional to the slope of the surface at that 
point. 

Consider first the sloping surface formed by a pair of oblique 
source— sink lines. The tip section lies along the lines of constant 
pressure of magnitude proportional to sinh~^ l/m. Between the 
sourco and sink the pressures are additive, so that 


q 


— m 
Jt ^ 1+m^ 


^ sinb~i ^ 
dx m 


w 


Ahead of or behind this section the pressvires cancel. 
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In case of a curved airfoil surface the chord can be divided 
into elements composed of source— sink pairs the strengths of which 
are proportional to the slope of the surface at the point in question. 
Each pair then contributes a pressure proportional to the local slope 
and contributes no pressure at other points. Hence, equation (l4) 
applies vrhen dz/dx is variable along the chord. 

Tine foregoing arguments of course apply only at the tip 
section of the oblique wing. At some distance from the tip section 
the overlapping isobars radiating from the tip again produce a 
quasi— cylindrical pressure field as in the case of the rectangular 
wing. Thus the resultant pressure distribution at either tip of a 
long oblique wing consists of two components, one given by equation 
(l4) and of the Ackeret type while the other component is equal to 
one-half the normal two-dimensional pressure distribution associated 
with the airfoil section. 


Figure 5 shows the pressures over a beveled-edge profile 
having 45"^ sweepback. The pressure distribution over the root 
section is given by 


^ = — m 
q It \Zi 4 ffl 2 


^ i Qo (x) - sinh ^ i Pq (x) 
dx (_ m 


(15) 


at a. great distance from either root or tip by 


^ = zi m 
q If 


^ Qo (x) 
dx 


(16) 


and at the tip by 


q 


It 


m 



Qo (x) + sinh ^ j Po (x) 


(17) 


To take account of the effect of compressibility we make use 
of the Prandtl transformation, increasing both the x dimensions 

1 

and the pressure coefficients by the factor \/ 1 - 142 , Replacing m 

by a/ 1-M- cot A, where a is the angle of sweepback, equation 
(l6) reduces to 
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^ 

p/2(V cos a)^ 


n 


1 

a/1— (M cos A)^ 


dz 

d(x cos a) 


Qo (x) 


(18) 


Thus at a great dista.nce from either root or tip the pressures follow 
a variation indicated hy the normal component of velocity V cos A. 

A.t the root section a component representing the Ackerot typo 
of pressvire distribution is added to equation (l8) . Ihis component 
is 


— 1 dz 

It V 1— (M cos a)'"^ d(x cos a) 


sinh ^ 


1 

^1—^12 cog ^ 


Po (x) 


(19) 


The factor sinlf^ ^ shows a logarithmic infinity at 

V l-f*I'' CCS A 

M = 1.0. Hence the pressure on the root section increases more 
rapidly with Mach, number than do the pressures at other sections 
of the swept—baclc wing. Fttrthermore, the sha.po of tho pressure 
distribution over the root section approaches the Ackeret shape 
more closely as tho Mach number approaches 1.0, As shown in 
reference 1 the pressure distribution on tho root section is exactly 
this shape at supei’scnic speeds , that is. 


Ap 

q cos2 A 


- 1 
^ \/l— (M cos a)^ 


dz 


d(x cos a) 


cosh~' 


VM*^— 1 cos Ay 


( 20 ) 


Since sinh~^ — >cosh~^ for large values of the argtanent, tho 
Bwept-back airfoil shows no discontinuity in the type of pressure 
distribution on passing througli the speed of sound. It will be 
evident that similar reasoning can be applied to the tip sections. 


AISFOIL OF BICONVEX SECTION 

Tho use of a. finite number of sources and sinks results in 
airfoil sections composed of straight segments. Such sections are 
undesirable, sinco thoy show infinite pressure peaks at the bonds 
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in tho surface. Surfaces having continuous curvature require 
continuous distribution of sources and sinks alined with the generators 
of the surface. The simplest of Hiese is the biconvex profile In 
wliich the upper and lower surfaces are parabolic arcs and have 
constant curvature. Such a. profile requires lino sources of finite 
strength to form the desired angles of intersection of the arcs at 
tho loading and trailing edges together with a unifonn distribution 
of sinks along the chord plane between the two sources. 

Tho pressure field for a uniform sheet of lino sources is 
obtained by integrating tho field of a single line source in tho x 
direction. This integral is 




sinh ^ 


^ dx = y sinh-^-^ - ^ y* 

|y»j m !yj m 


sinh ^ 



( 21 ) 


The integration for a source sheet is actually somewhat simpler if 
the interference of a bilaterally symmetrical arrangement of sources 
is considered simifltaneously . Tho influence of the symmetrical, or 
conjugate, arrangement is obtained by substituting -m for m in 
equation (21) . Denoting x-n:ny by and y+mx by y* we have 


F ' f 


^sinh' 


— 1 


+ sinh' 


— 1 



1 

ra 



sinh ^ 



sinh ^ 



(22) 


To obtain a complete swept— back wing it is necessary to add 
a number of component pressure fields as oxpla.inod in reference 1. 
For an infinite swoptMsack wing with leading and trailing edges at 
yf = +111 and -m, I'espectively, on one side, and at y' = +m 
and -m, reBpoctivol5’, on the other side, there is obtained 


4^ _ 2 B1 

1 ^ l-Hn i^ 


' I 2LL {' sinli ^ 


x'+l 

iy*-m| 


sinh ^ 


+ 


Zl 

m 


sinli ^ 


x’-l 

!y'-^ 


sinh ^ 



+ 2 




+ 2 Qi 




(23) 
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vhere thickness— chord ratio of the biconvex profile. 

The terms Qi 
biconvex airfi 

terms is the result of the assumption that the tips are n.'emoved to 
a great distance . 


(• 


ZL) 
m ) 


represent the pressiaro distribution on the 


3il In two-dimonsional flow. The appearance of these 


At the root section (y = 0) equation (23) reduces to 


Ap _ 2 ra 
q Jt iJl+TSl 2 


eXax 


4 Qi (x) - 4 sinli~^ 1 Pi (x) 


(24) 


Figure 6 shows pressure distributions at various stations along 
the span foi’ a biconvex wing with 60® swoopback. The curves assume 
the two-dimensional form at a relatively short distance (y - c) 
from the root section, and similar behavior is to bo expected near 
the tips. Hence the assumirtion of infinite aspect ratio should 
apply very nearly at any section situated more than one— half chord 
length from either root or tip. 


Figiuro T shows the effect of I-Iach number on the pressures 
over the root section and illustrates the progressive change to the 
supersonic type as the Mach nvimber approaches 1,0. It can be seen 
that an increase in Mach number will not only increase the dlstoio- 
tion of the pressure distribution but will increase the extent of 
the distortion along the span. 


An interesting point to be noted is that not all sections of 
the swept-back wing have zero pressure drag. A positive drag appears 
on the root sections and a negative drag on the tip sections. Hence 
the spanwiso drag distribution is qualitatively similar to that at 
supersonic speeds though, of course, the not subsonic pressure dra.g 
is zero. 


Ames Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Moffett Field, Calif., May 194?. 
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APPENDIX 

Sj^bols 

V flight velocity 

M Mach number 

x,y,z coordinates 

I point on x— axis 

Ti point on y-axis 

cp disturbance, velocity potential 

u,v,-w disturbance velocity components 
p local pressure 

q dynamic pressure (^pVS) 

p air density 

Pn,Qn Legendre functions 

D differential operator (d/dx) 

t thickness of wing 

c chord of wing (measured along x) 

m slope of line soiorce (absolute value) 

X* X + my 

y* y - mx 

X* X — my 

y' y + mx 

E.P. Beal part 
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Fig. 3 



F/CU^^E 3 C0AJ3-TEUCT10E' OF THE PEE3SUR.E 


D/STP/Sur/OAJ Ol/EE A BEVELED LEADIK>C EDGE. 


Fig. 4 
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Fig. 7 
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